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Abstract 

Stability of droplets of ud matter, viz udilets, is studied using non- 
relativistic and relativistic bag models. It is found that there is a 
strong possibility for forming hypercharged metastable udilets with 
charge fraction Z/A > 1 in semicentral relativistic heavy ion collisions. 
It is then speculated that protons hadronised from these hypercharged 
udilets may be contributing towards the directed sideward flow ob- 
served recently in experiment. This exotic state with large charge to 
mass ratio presents an unique possibility for the spontaneous emission 
of positrons due to the 'sparking' of the vacuum. 
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The principle goal of present day relativistic heavy ion collision experi- 
ments is to find new states of matter such as quark gluon plasma Such 
pursuits are based on the viewpoint that there exists a transition between 
nuclear and quark matter [2-12]. Logically, one would expect the transition 
to take place when the nuclear matter is compressed to a density at which 
nucleons overlap to a significant extent. This density is roughly given by 
Pn ~ 3/ (47rr]y), where is the nucleon size. Taking ~ 0.8 — 1.0 fm, we 
get pn ~ 0.24—0.40 fm~ 3 . Consequently, it has been conjectured ||, [7|, [TT], [12 



that a deconfined quark phase(DQP) could be produced in heavy-ion colli- 
sions(HIC). 

Over the years, considerable amount of work has been done to find 'un- 
ambiguous' signatures of DQP]13[. In this regard, probable formation and 
detection of small lumps of strange quark matter in relativistic HIC is being 
promoted as an unambiguous signature [[l^]. This speculation is based on 
the fact that in ideal baryon rich quark matter in thermodynamical equilib- 
rium , the production of ss pairs is enhanced compared to the production 
of light quarks (u & d) as soon as the Fermi energy of the already present u 
and d quarks is higher than twice the strange quark mass. If the tempera- 
ture and density attained in a semi central/peripheral collision is such that 
the production of antiparticles and strange quarks are suppressed, then the 
chances are that a blob of DQP consisting of only u and d quarks, namely 
an udilet, will be formed. The probability of forming an udilet decreases 
as the collision becomes either very central or very peripheral. It would be 
interesting to understand how does such a state show itself, and if there is 
any specific signal, quite different from those characteristic of very central 



collisions. In this work, we mainly focus upon these aspects and investigate 
them using a non-relativistic model(NRM) [|14j] as well as the usual relativistic 
bag model(RBM)p|. 

Recently, we have developed |14| a consistent mass formula for multiquark 
droplets based on a non-relativistic framework. This formula is quite success- 
ful in describing the properties of strangelets, and also is found to be fully con- 
sistent with relativistic studies. Omitting the details, we give here only the 
final expression. Consider a system of N u and number of u and d quarks 
respectively, with mass m, confined in a spherical box of size R. For a given 
baryon number A, we define an asymmetry parameter, 5 = (N^ — N u )/N q , 
where N q = N u + = 3A. Then the quark flavor densities are defined as 
p u = 0.5p O(J (l — 5) and pd = 0.5p og (l + 5), where p oq is the equilibrium total 
quark number density. We can then calculate, for a given 5, the total mass 



energy of an udilet using the expression [[14 
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given in terms of the equilibrium baryon number density p b(=Poq/3) and the 
total charge Z = 0.5 A(l — 36). Here, we take the constituent quark mass m 
to be 300 MeV|]l4|- We also give the corresponding expression obtained in 
the case of the usual relativistic bag model where the u and d quarks are 
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massless. 



M RL v(A) = 




It may be noted that unlike in the non-relativistic model, the surface energy 
coefficient for an udilet is zero in the present case. Here, we have omitted the 
lowest order quark-quark interaction with the supposition that it can effec- 
tively be absorbed into the bag parameter, or equivalently into the saturation 
density p Q b as has been done by Berger and Jaffe||. 

We then calculate the mass per baryon of an udilet at the equilibrium 
density p Q b for various values of charge fraction fz = Z/A. We have taken 
p b = 0.25 fm~ 3 in the case of NRM. And in the case of RBM, the value 
of p b is taken to be 0.40/m~ 3 so that M RL y(p ob )/A > M N , where Mtv=939 
MeV. Results so-obtained using Eqs.(l) and (2) are shown in Fig.l for a 
representative value of A. It is seen that the mass per baryon shows a well 
defined minimum at /™ n ~ 0.4 for A = 300 in both the models. Further, 
we find in both the models, f™ m increases as A is decreased and/or p ob is 
increased. It is satisfying to find that our non-relativistic model qualitatively 
agrees with the relativistic bag model. 

Let us now suppose that an udilet with fz — 0.4 and A ~ 300 is formed 
in a semi central/peripheral HIC at relativistic energies. The two important 
channels for the decay of the udilet are fission and strong nucleon emission. 
Firstly, let us investigate whether the udilet undergoes fission. For fission to 
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be energetically favourable, we must have the mass energy difference Ef — 
M(A, Z) — Mi(Ai, Zi) — M 2 (A 2 , Z 2 ) > 0. Considering symmetric fission and 



a liquid drop model like expansion for the mass energy [Eq. (13) of Ref.fl4 
one can rewrite the criterion Ef > as Z 2 /A > 0.26a s /0.37ac where a s and 
ac are respectively the surface and Coulomb energy coefficients. In the case 
of RBM, where a s ~ 0, the criterion Ef > can be expressed in terms of the 
curvature energy coefficient 

Z 2/ A 2/3 > o.59a cv /0.37a c . Further, we are 
aware from the nuclear fission studies that eventhough Ef > 0, instantaneous 
fission does not always happen. An additional criterion for such a decay 
process is that Ef/Eb > 1, where Ej, = ZiZ 2 e 2 /[r (Ay 3 + A 2 ^ 3 )] is the 
Coulomb barrier energy and r = (3/(47rp {,)) 1//3 . In Fig. (2), we have plotted 
Ef/Eb calculated using Eqs.(l) and (2) as a function of Z/A, for A = 300 
and A = 40. In the NRM, for A = 40, fission is not allowed. For A = 300, 
Ef/Eb being negative upto fz — 1.6, fission is disallowed and over the range 
fz ~ 1.6 — 2.0, though Ef/E b is slightly positive, instantaneous fission is 
disallowed. But, in the RBM, udilets with large value of A are susceptible 
towards fission, which may be due to the vanishing of a s . For example, 
Ef/Eb ~ 1 at f z ~ 0.7, for A = 300; and at f z ~ 1.0, for A = 200. Further, 
for a given fz, the ratio Ef/Eb becomes more negative as A decreases and/or 
p b increases. However, in the case of relativistic bag model as the finite-size 
coefficients such as curvature a cv and Coulomb a c are proportional to p^ 3 , 
the ratio Ef/Eb is independent of p Q b- Since, an udilet formed in a HIC is 
most likely to have A < 300 and fz — 0.3 — 0.5 and fission is disallowed for 
these values of A and fz in both the models, the predominant channel for 
decay is the strong nucleon emission. 



Among the baryons like neutron, proton and A's, the neutron emission is 
the most favourable one because of low mass and non-existence of Coulomb 
barrier. As an udilet with fz ~ 0.4 and A ~ 300 emits neutrons, it is driven 
towards higher values of fz- This process is energetically allowed as long as 
Q n = M(A, Z) - M(A - 1, Z) - M N > 0. We find that once the charge 
fraction fz reaches a value of about 1.2 with A ~ 102 and 1.9 with A ~ 64 
in the cases of RBM and NRM respectively, the neutron emission stops. It 
is also seen that the limiting values of fz increases as A decreases and/or 
Pol, increases. It maybe noted that the limiting value of fz obtained in the 
case of RBM is lower than the one obtained in the case of NRM. At this 
stage of hyper- char geness, the system is unstable against proton, A + and 
A ++ decays. Here, the proton emission is the most favourable one, and is 
allowed as long as the corresponding Q— value, Q p > 0. We find in the case 
of NRM, the proton emission continues upto fz ~ 2 with A ~ 58, where the 
entire udilet is composed of only u quarks. This non-relativistic system is 
highly unstable against strong A ++ decay, and thus it finally disintegrates. 

In the case of RBM, the proton emission continues upto fz ~ 1-9 with 
A ~ 20. At this stage, the system is unstable against strong A ++ decay. 
However, it must be noted that for such small baryon numbers, shell effects 



are important [[16] and shall decide the stability of the udilet. For this reason, 
we recalculated the Q- value, Q A ++ = M(A, Z) - M(A - 1, Z - 2) - M A+ +, 
using the shell model calculations [16|, [H]. The total mass energy M(A,Z) 
of an udilet at the saturation point is obtained using the expression, 
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M(A,Z) sR 



N u N d 

Qi^i + 9iVi + 0.6Z 2 e 2 

i=l i=l 



(3) 



where R = tqA 1 ^ 3 and gi is the degeneracy factor. It may be noted that the 
values of A and f z of the udilet remaining after the emission of neutrons 
and protons depend upon the characteristics^, f z , p b) acquired at its birth. 
For example, if the initial value of A is 200 instead of 300, then the limiting 
values at the end of neutron and proton emissions are A ~ 14 and fz ~ 2. 
Therefore, it may be appropriate to study the dependence of Qa++ on A, fz 
and p ob . 

To do so, we calculated Qa++ and the Coulomb barrier energy Vcb = 
2{Z — 2)e 2 (A — l) _1 ^ 3 /r for various values of A, keeping fixed p ob and f z . 
Results are displayed in Figs. (3) and (4) for two values of fz and three 
values of p ob . It is interesting to see that certain udilets are absolutely stable 
against A ++ decay, i.e. Qa++ < 0. For a given A, the stability increases 
as fz and p Q b decreases. Under these conditions, udilets can decay only by 
converting some of the u quarks into d quarks through weak leptonic decay 
process. Hence, this study predicts quite a good possibility for the formation 
of metastable udilets of relatively long decay lifetime, which is likely to be 
observed in experiments. The characteristic signature of such udilets is their 
large charge to mass ratio, Z/A > 1, which is in sharp contrast to the low 
value Z/A « 1 found in case of strangelets. 

Thus, from both the relativistic and non-relativistic studies, an optimistic 
scenario about DQP in form of udilets with large charge to mass ratio as char- 
acteristic signature emerges. While in former case metastable udilets with 
relatively large lifetime are formed, in the latter case they disintegrate by 
proton, A + and A ++ decay processes. The predominance of proton in this 
phenomenon is quite evident. This is suggestive of the interesting possibility 
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that the 'direct' protons as well as those from A ++ — > p + tt + decays, orig- 
inating from the udilets, may be contributing toward the directed sidewards 
flow observed recently [18|]. In this context, the observations that protons 
are the major carriers of flow and the magnitude of the flow decreases as 
the collision becomes very central are quite inspiring to us. Further, in the 
present study we have shown that the udilets formed after neutron emis- 
sions has fz — 1-2 — 1.9, depending upon the model. This presents a novel 
possibility of a system where charge number substantially exceeds the mass 
number. This then provides a highly favourable situation for reaching the 
critical value for Z (~ 137) for the occurence of the spontaneous emission of 
positrons due to the 'sparking' of the vacuum |TI| . Finally, we would like to 
comment that for the aforesaid possibilities to be realised, the initial param- 
eters of an udilet, i.e. A, fz and p Q b, should be such that the charge fraction 
acheived at the end of neutron emissions is greater than one. 

In conclusion, stability of udilets is studied using non-relativistic and rel- 
ativistic bag models. It is found that there is a strong possibility for forming 
hypercharged metastable udilets with charge fraction Z/A > 1 in semicentral 
relativistic heavy ion collisions. This large value of charge fraction can serve 
as an unique signature for DQP. This is in sharp contrast to the strangelets 
which are characterised by small charge fraction. It is then speculated that 
protons hadronised from these hypercharged udilets may be contributing 
towards the directed sideward flow observed recently in experiments. An ex- 
citing novel possibility of forming a system having large charge number with 
relatively small baryon number raises new hope of the spontaneous emission 
of positrons due to the 'sparking' of the vacuum. 
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FIGURE CAPTIONS 



FIG.l : Total mass per baryon of an udilet calculated using Eqs.(l) [NRM] 
and (2)[RBM] are shown as a function of the charge fraction Z/A for baryon 
number A = 300. 

FIG. 2 : Ratio of the fission energy Ef to the Coulomb barrier energy 
calculated using Eqs.(l)[NRM] and (2)[RBM] are shown as function of the 
charge fraction Z/A for two values of baryon number A. 

FIG. 3 : Values of /Vcb obtained from shell model calculations[Eq.(3)] 
are shown as a function of baryon number A for three values of saturation 
baryon number density p &, keeping fixed the charge fraction fz = Z/A at 
1.8. 

FIG.4 : Same as Fig.3, but for f z = 1.3. 
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